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Abstract 

We discuss the critical behavior of several three-dimensional mag- 
netic systems, such as pure and randomly dilute (anti)ferromagnets 
and stacked triangular antiferromagnets. We also discuss the nature 
of the multicritical points that arise in the presence of two distinct 
0(n)-symmetric order parameters and, in particular, the nature of 
the multicritical point in the phase diagram of high-T c superconduc- 
tors that has been predicted by the SO (5) theory. For each system, 
we consider the corresponding Landau-Ginzburg- Wilson field theory 
and review the field-theoretical results obtained from the analysis of 
high-order perturbative series in the frameworks of the e and of the 
fixed-dimension d = 3 expansions. 



1 Introduction 



In the framework of the renormalization-group (RG) approach to critical 
phenomena, a quantitative description of many continuous phase transitions 
can be obtained by considering an effective Landau- Ginzburg- Wilson (LGW) 
theory, containing up to fourth-order powers of the field components. The 
simplest example is the 0(iV)-symmetric <p A theory, 



n 



O(N) 



J d d x [± J2(d^) 2 + \r £ *? + i« £ S?*?] , (1) 

i i ' ij 



where $ is an iV-component real field. This model describes several univer- 
sality classes: the Ising one for N = 1 (relevant for the liquid-vapor transi- 
tion in simple fluids, for the Curie transition in uniaxial magnetic systems, 
etc.), the XY one for N = 2 (relevant for the superfluid transition in 4 He, 
and the magnetic transition in magnets with easy-plane anisotropy, etc.), 
the Heisenberg one for N — 3 (it describes the Curie transition in isotropic 
magnets), the 0(4) universality class for N = 4 (relevant for the finite- 
temperature transition in two-flavor quantum chromo dynamics, the theory 
of strong interactions); in the limit N — > it describes the behavior of dilute 
homopolymers in a good solvent in the limit of large polymerization. See, 
e.g., Refs. PI |2] for recent reviews. But there are also other physically in- 
teresting transitions described by LGW field theories characterized by more 
complex symmetries. 

The general LGW Hamiltonian for an ^-component field $j can be writ- 
ten as 

n = /^[^ECW' + iE^ + ^E^^i^]. ( 2 ) 

i i ijkl 

where the number of independent parameters and u^ki depends on the 
symmetry group of the theory. An interesting class of models are those in 
which J2i ®i is the only quadratic invariant polynomial. In this case, all T{ 
are equal, = r, and Uiju satisfies the trace condition jSj 

J2u iiM (x5 k i. (3) 

i 

In these models, criticality is driven by tuning the single parameter r. There- 
fore, they describe critical phenomena characterized by one (parity-symmetric) 
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relevant parameter, which often corresponds to the temperature. Of course, 
there is also (at least one) parity-odd relevant parameter that corresponds 
to a term YU that can be added to the Hamiltonian (J2J). For symmetry 
reasons, criticality occurs for hi — *■ 0. There are several physical systems 
whose critical behavior can be described by this type of LGW Hamiltoni- 
ans with two or more quartic couplings, see, e.g., Refs. |UI2]. More general 
LGW Hamiltonians, that allow for the presence of independent quadratic 
parameters 7*j, must be considered to describe multicritical behaviors arising 
from the competition of distinct types of ordering. A multicritical point can 
be observed at the intersection of two critical lines with different order pa- 
rameters. In this case the multicritical behavior is achieved by tuning two 
relevant scaling fields, which may correspond to the temperature and to an 
anisotropy parameter [5]. 

In the field-theoretical (FT) approach the RG flow is determined by a 
set of RG equations for the correlation functions of the order parameter. In 
the case of a continuous transition, the critical behavior is determined by the 
stable fixed point (FP) of the theory, which characterizes a universality class. 
The absence of a stable FP may instead be considered as an indication for a 
first-order transition, even in those cases in which the mean-field approxima- 
tion predicts a continuous transition. But, even in the presence of a stable 
FP, a first-order transition may still occur for systems that are outside its 
attraction domain. The RG flow can be studied by perturbative methods, 
by performing an expansion in powers ofe = 4 — djUJora fixed- dimension 
(FD) expansion in powers of appropriate zero-momentum quartic couplings 
[7j. It is also possible to combine the FD approach with the 4 — d regular- 
ization, obtaining the minimal-subtraction scheme without e expansion [S]. 
In these perturbative approaches, the computation and the resummation of 
high-order series is essential to obtain reliable results for three-dimensional 
transitions (see Refs. PU 12] for reviews). Beside improving the accuracy, in 
some cases high-order calculations turn out to be necessary to determine the 
correct physical picture in three dimensions. 

In this paper we discuss the applications of the FT approach to the study 
of several three-dimensional transitions in magnetic systems. Beside the 
0(iV)-symmetric Hamiltonian 7io(N) that is relevant for magnetic transition 
in uniaxial, easy-plane, and isotropic magnets, we consider LGW Hamil- 
tonians characterized by more complex symmetries: in particular, we con- 
sider Hamiltonians invariant under the cubic group, under 0(M)<S)0(N) and 
0(M)©0(iV). They allow us to investigate the effects of the addition of un- 
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correlated nonmagnetic impurities to magnetic systems, the critical behavior 
of frustrated spin models with noncollinear order — we can check whether a 
new chiral universality class exists — and the nature of the multicritical points 
arising from the competition of distinct 0(n) order parameters. We present 
an overview of results obtained by using perturbative FT methods and com- 
pare them with other theoretical approaches that provide results of compa- 
rable precision — first of all, with results obtained using lattice techniques, 
e.g., Monte Carlo and high-temperature expansions — and with experiments. 

The paper is organized as follows. Sec. El is dedicated to the standard 
O(N) universality classes. We present a brief overview of FT results for 
the universality classes that are relevant for magnetic systems, such as the 
Ising, XY, and Heisenberg ones. In Sec. |3] we analyze the stability of the 
0(N) -symmetric FP under generic perturbations in three-dimensional N- 
component systems. This analysis is useful to assess the relevance of magnetic 
anisotropies in the critical behavior of physical systems. In Sec. 01 we discuss 
the effects of the addition of quenched nonmagnetic impurities on the critical 
behavior of magnetic systems, focusing on pure systems that have transitions 
belonging to the standard 0(N) universality classes. In Sec. E]we consider 
frustrated spin models with noncollinear order, whose critical behavior is 
effectively described by an 0(M)®0(jV)-symmetric Hamiltonian with M = 
2. These models are relevant for stacked triangular antiferromagnets. In 
Sec. Elwe consider the 0(n 1 )©0(n 2 )-symmetric 4 theory that describes the 
universal behavior near the multicritical point where two critical lines with 
symmetry O(ni) and 0(7x2) meet. Finally, we discuss the relevance of these 
results for anisotropic antiferromagnets in the presence of a magnetic field 
and for high-T c superconductors within the so-called SO (5) theory. 

2 The O(N) universality classes 

The O(N) universality classes are relevant for magnetic systems, because 
they describe the critical behaviors in uniaxial, easy-plane, and isotropic 
magnets, corresponding respectively to the Ising (N = 1), XY (N = 2) 
and Heisenberg (N = 3) universality classes. The universal critical proper- 
ties can be investigated by using perturbative FT methods applied to the 
0(iV)-symmetric Hamiltonian (£[]). The oldest perturbative method is the e 
expansion in which the expansion parameter is e = 4 — d [6 . Subsequently, 
Parisi [7] pointed out the possibility of using pertubation theory directly at 
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the physical dimensions d = 3 and d = 2 in the massive (high-temperature) 
phase. 



2.1 The fixed-dimension expansion 

In the FD expansion one works directly in d = 3 or d = 2. In this case 
the theory is super-renormalizable since the number of primitively divergent 
diagrams is finite. One may regularize the corresponding integrals by keeping 
d arbitrary and performing an expansion ine = 3 — <i or e = 2 — d. Poles in e 
appear in divergent diagrams. Such divergences are related to the necessity 
of performing an infinite renormalization of the parameter r appearing in the 
bare Hamiltonian, see, e.g., the discussion in Ref. This problem can be 
avoided by replacing the quadratic parameter r of the Hamiltonian with the 
mass m (inverse second-moment correlation length) defined by 



o i 0r< 2 V 

m 



r(2)(o) dp 2 



(4) 

p 2 =0 



where the function T^^p 2 ) is related to the one-particle irreducible two-point 
function (i.e. the inverse two-point function of the order-parameter field) by 

r2 } (p) = ^r (2 V)- (5) 

Perturbation theory in terms of m and u is finite in d < 4. The critical limit 
is obtained for m — ► 0. To handle it, one considers appropriate RG functions. 
Specifically, one defines the zero-momentum four-point coupling g and the 
field-renormalization constant by 

T^(p) = 6 ab Z^[m 2 +p 2 + 0(p')}, (6) 
ritL(O) = Z^m A - d9 - (5 ab 5 cd + S ac 5 M + 5 ad 5 bc ) , (7) 

where an are one-particle irreducible correlation functions. Then, one 
defines a coupling-renormalization constant Z g and a mass-renormalization 
constant Z t by 

u = m 4 ~ d gZ g Z^ } ril' 2) (0) = 5 ab Z~\ (8) 

where r„ b ' (p) is the one-particle irreducible two-point function with an in- 
sertion of |0 2 . The renormalization constants are determined as perturbative 
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expansions in powers of g. The FP of the model is determined by the non- 
trivial zero g* of the /^-function 



(d - A)g 



i + ^MW) 



Then, one defines 



Vt{g) ■ 



dhiZ, 



dlnm 
d\nZ t 



d\nm 



--M 



d\nZ 4 

dg 
d\nZ t 



dg 



Finally, the critical exponents are given by 



V 
v 

00 



[2 - V ^g*) + Vtigir 1 



(9) 

(10) 
(11) 



(12) 
(13) 
(14) 



where uo is the exponent associated with the leading irrelevant operator. All 
other exponents can be obtained using the scaling and hyperscaling relations. 

In three dimensions the /3-function is known to six loops ^U], while the 
RG functions associated with the critical exponents have been computed to 
seven loops [TT] . FT perturbative expansions are divergent. Thus, in order 
to obtain accurate results, an appropriate resummation is required. This 
can be done by exploiting their Borel summability, and the knowledge of the 
high-order behavior of the expansion, see, e.g., Ref. in Table^we report 
some recent results for the critical exponents of the Ising, XY, and Heisen- 
berg universality classes. For comparison, we also report estimates from lat- 
tice techniques, such as high-temperature (HT) expansions and Monte Carlo 
(MC) simulations. The agreement is very good. A much more complete list 
of results and references can be found in Ref. Of course, it is important 
to note that there is also a good general agreement with experiments, see, 
e.g., the list of results reported in Ref. 0. 

Since this FD expansion is based on zero-momentum renormalization con- 
ditions, it is not suited for the study of the low-temperature (broken) phase. 
In order to investigate the low-temperature phase other perturbative FD 
schemes should be considered, such as the minimal-subtraction scheme with- 
out e expansion [S]. 
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Table 1: Three-dimensional estimates of the critical exponents from the anal- 
ysis of the fixed-dimension (FD) and e expansions, and from lattice techniques 
such as high-temperature (HT) expansions and Monte Carlo simulations. 



N 


1/ 


V 


a 


7 





1 FD exp [T3] 
e exp 22] 
HT exp [H] 

mc 031 


0.6303(8) 
0.6305(25) 
0.63012(16) 
0.6297(5) 


0.0335(25) 
0.0365(50) 
0.03639(15) 
0.0362(8) 


0.1091(24) 
0.108(7) 
0.1096(5) 
0.1109(15) 


1.2396(13) 
1.2380(50) 
1.2373(2) 
1.2366(15) 


0.3258(14) 
0.3265(15) 
0.32653(10) 
0.3262(4) 


2 FD exp [HI 
e exp 22] 
HT exp [El 

MC 16 


0.6703(15) 
0.6680(35) 
0.67155(27) 
0.6716(5) 


0.0354(25) 
0.0380(50) 
0.0380(4) 
0.0380(5) 


-0.011(4) 
-0.004(11) 
-0.0146(8) 
-0.0148(15) 


1.3169(20) 
1.3110(70) 
1.3177(5) 
1.3177(10) 


0.3470(11) 
0.3467(25) 
0.3485(2) 
0.3486(3) 


3 FD exp [T21 
e exp ^21 
HT exp [32J 

MC P2J 


0.7073(35) 
0.7045(55) 
0.7112(5) 
0.7113(11) 


0.0355(25) 
0.0375(45) 
0.0375(5) 
0.0378(6) 


-0.122(10) 
-0.114(16) 
-0.1336(15) 
-0.1339(33) 


1.3895(50) 
1.382(9) 
1.3960(9) 
1.3957(22) 


0.3662(25) 
0.3655(50) 
0.3689(3) 
0.3691(6) 



2.2 The e expansion 

The e expansion jB] is based on the observation that, for d = 4, the theory 
is essentially Gaussian. One considers the standard perturbative expansion, 
and then transforms it into an expansion in powers of e = 4 — d. In practice, 
the method works as in the FD expansion. One first determines the expansion 
of the renormalization constants Z g , Z^, and Z t in powers of the coupling g. 
Initially, they were obtained by requiring the normalization conditions ©, 
(0), and (|5ji . However, in this framework it is simpler to use the minimal- 
subtraction scheme [T%j . Once the renormalization constants are determined, 
one computes the RG functions (3(g), r]$(g), and rj t (g) as in Sec. 12.11 The 
FP value g* is obtained by solving the equation /3(g*) = perturbatively in 
e. Once the expansion of g* is available, one obtains the expansion of the 
exponents, by expressing r]^(g*) and T] t (g*) in powers of e. In this scheme, 
five-loop series for the exponents were computed in Refs. [121 [20]. The results 
of their analysis, exploiting Borel summability, are reported in Tabled Again 
the results are in good agreement with the other approaches. 
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3 Stability of the O(N) -symmetric fixed point 
in 3-D multicomponent systems 

In order to discuss the stability of the 0(N) FP in a generic iV-component 
system, we consider the general problem of an 0(iV)-symmetric Hamiltonian 
in the presence of a perturbation P, i.e., 



n 



J d d x [\{d^f + ±r$ 2 + i M ($ 2 ) 2 + h p P] , (15) 



where $ is an iV-component field and h p an external field coupled to P. 
Assuming P to be an eigenoperator of the RG transformations, the singular 
part of the Gibbs free energy for the reduced temperature t — > and h p — ► 
can be written as 

F^(t,h J ,)a\t\*'?(h J ,\t\-*>), (16) 

where (fi p = y p v is the crossover exponent associated with the perturbation 
P, y p is the RG dimension of P, and J-(x) is a scaling function. If y p > the 
pertubation is relevant and its presence causes a crossover to another critical 
behavior or to a first-order transition. 

In order to discuss the stability of the 0(N) FP in general, we must con- 
sider any perturbation of the 0(N) FP. We shall first consider perturbations 
that are polynomials of the field $ a . Any such perturbation can be written 
|2*T] as a sum of terms P^ 1 i'"' ,ai , k > I, which are homogeneous in $ a of degree 
k and transform as the Z-spin representation of the 0(N) group. Explicitly, 
we have 

p ai,..., ai = ^2j(*-0/2g«i,..,«, 

where Q«'-' tt « is a homogeneous polynomial of degree I that is symmetric 
and traceless in the I indices. The lowest-order polynomials are 

Qi = $ a , (18) 

Qab = $ a $ fe _ 1 S ab$2 ^ 

qabc = $ a $ 6 $c _ j^-^> 2 (<S> a 5 bc + <5> b 5 ac + $ c <5 ab ) , (20) 

qabcd _ ^a^b^c^d 
1 



N + 4 

+ __^___ (<|)2)2 + 5 ac 5 bd + ^60) (21) 



7 



The classification in terms of spin values is particularly convenient, since 
polynomials with different spin do not mix under RG transformations. On 
the other hand, operators with different k but with the same I do mix under 
renormalization. At least near four dimensions, we can use standard power 
counting to verify that the perturbation with indices m, I mixes only with 
Pk',h k' < k. In particular, P^ renormalizes multiplicatively and is therefore a 
RG eigenoperator. Moreover, if y k j is the RG dimension of the appropriately 
subtracted P^, one can verify that for small e, y^i < 0, for I > 5, i.e. the 
only relevant operators have I < 4. It is reasonable to assume that this 
property holds up to e = 1. In principle, one should also consider terms with 
derivatives of the field, but again one can show that they are all irrelevant or 
redundant. Note that the condition that &l is the only quadratic invariant 
of the theory forbids the presence in the Hamiltonian of any spin-2 term Pf 2 . 
Analogously, the trace condition Q forbids quartic polynomials transforming 
as the spin-2 representation of the 0(N) group, i.e. the operators P£ b 2 . 

Let us first discuss the case of perturbations that are even under parity. 
Beside the 0(iV)-symmetric terms P2 i o = *£ 2 an d P^o — ($ 2 ) 2 there are 
only three other polynomial perturbations that must be considered, i.e., P 2 b 2 , 
P42, and Pf 3 ^} According to the above-reported general arguments, the 
stability properties of the 0(N) FP under even perturbations can be obtained 
by determining the RG dimensions y^j of the five classes of operators P2,o, 
P 22 , Pifi, P£ b 2 and P^ d . Of course, the dimensions y^ do not depend on 
the specific values of the indices and thus one can consider any particular 
combination for each class of operators. In Table 121 we report FT estimates 
of the RG dimensions yk,i for N = 2, 3, 4, 5, obtained from the analysis of six- 
loop FD and five-loop e series |2*2*1 1231 IT3*] . 2 The RG dimensions of P2,o and 
of Pa. can be computed directly in the O(N) -invariant theory. In particular, 
1/2,0 = 1/v and 2/4,0 = —00, where u is the leading irrelevant exponent in the 
0(N) -invariant theory. The quadratic perturbations P 2 a 2 are relevant for the 
description of the breaking of the O(N) symmetry down to 0(M)©0(iV— M). 
Since y 2i2 > 0, they are always relevant. The RG dimension y^ 2 is negative for 

Note that, according to the above-reported discussion, Pf 2 and P£\ cd are RG eigen- 
operators, while P4 \ must be in general properly subtracted, i.e. the RG eigenoperator 
is P4 2 + zP§2 f° r a suitable value of z. The determination of the mixing coefficient z 
represents a subtle point in the FD expansion [7j , but is trivial in the MS scheme in 4 — e 
dimensions, in which operators with different dimensions never mix so that z — 0. 

2 Results obtained in other theoretical approaches and in experiments can be found in 
Rcfs. and references therein. 
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Table 2: Three-dimensional estimates of the RG dimensions y^j from e and 
FD expansions, from Ref. [T3] for ?/ 2 ,o an d 2/4,0, from Ref. [22] for 2/2,2, from 
Ref. [231211 for 1/4,2 and 2/ 4j4 , and from Ref. [22] for 7y 3)3 . 



TV 




2/2,0 = 


2/2,2 


2/4,0 


2/4,2 


2/4,4 


2/3,3 


2 


FD 

e 


1.493(3) 
1.497(8) 


1.766(18) 
1.766(6) 


-0.789(11) 
-0.802(18) 


-0.624(10) 


-0.103(8) 
-0.114(4) 


0.897(15) 
0.90(2) 


3 


FD 

e 


1.414(7) 
1.419(11) 


1.80(3) 
1.790(3) 


-0.782(13) 
-0.794(18) 


-0.550(14) 


0.013(6) 
0.003(4) 


0.96(3) 
0.97(4) 


4 


FD 

e 


1.350(11) 
1.357(15) 


1.82(5) 
1.813(6) 


-0.774(20) 
-0.795(30) 


-0.493(14) 


0.111(4) 
0.105(6) 


1.03(3) 
1.04(5) 


5 


FD 

e 


1.312(12) 
1.333(36) 


1.83(5) 
1.832(8) 


-0.790(15) 
-0.783(26) 


-0.441(13) 


0.189(10) 
0.198(11) 


1.07(2) 
1.08(4) 



any N, so that the corresponding spin-2 perturbation P 4 \ is always irrelevant. 
On the other hand, the sign of 7/4,4 depends on N: it is negative for N = 2 
and positive for N > A. For iV = 3 it is marginally positive, suggesting the 
instability of the 0(3) FP under generic spin-4 quartic perturbations. 

Let us now briefly consider the odd perturbations. The RG dimension 
yi t i of is given by y 1; i = (f3 + y)/v = 5/2 — ^/2. The operator P^ 
is redundant, because a Hamiltonian term containing P^i can be always 
eliminated by a redefinition of the field $ a . The RG dimension 7/3 3 of the 
spin-3 perturbations P3 | c has been computed to six loops in the FD expansion 
and to five loops in the e expansion [23]. The analyses of the series provide 
the estimates reported in Table [21 which show that the spin-3 perturbation 
P33 is relevant for any N. 

The stability of the O(N) FP under spin-4 perturbations can also be 
inferred from the RG flow of the cubic-symmetric LGW Hamiltonian for an 
iV-component field 

r ( 1 N 1 N N } 

*. = J*° {IX [( W + r*!\ +^[ w (E $ i) 2 + v E $ i]|- ( 22 ) 

The cubic-symmetric LGW Hamiltonian 7i c is physically relevant because 
it takes into account the most important source of magnetic anisotropy in 
Heisenberg systems. Indeed, the magnetic interactions in crystalline solids 
with cubic symmetry like iron or nickel are usually modeled by using the 
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0(3)-symmetric Heisenberg Hamiltonian. However, this is a simplified model, 
since other interactions are present. Among them, the magnetic anisotropy 
that is induced by the lattice structure (the so-called crystal field) is particu- 
larly relevant experimentally, see, e.g., Ref. j2H]. I n cubic-symmetric lattices 
it gives rise to additional single-ion contributions, the simplest one being 
X)i s*j 4 - These terms are usually not considered when the critical behavior of 
cubic magnets is discussed. However, this is strictly justified only if these 
nonrotationally invariant interactions, that have the reduced symmetry of 
the lattice, are irrelevant in the RG sense. This issue can be investigated by 
considering the cubic-symmetric Hamiltonian TC C |4 . Its RG flow has been 
much studied using various FT methods [2]. The O(N) FP turns out to be 
unstable for N > N c with iV c 3. The most accurate results have been 
provided by analyses of high-order FT perturbative expansions, six-loop FD 
and five-loop e series, see, e.g., Refs. j2H 123 12H] : they find N c < 2.9 in three 
dimensions and the existence of a stable FP characterized by a reduced cubic 
symmetry for iV > N c . Therefore, in three-component magnets the isotropic 
FP is unstable, and the RG trajectories flow toward a stable cubic FP, see, 
e.g., Ref. |21]. However, differences between the Heisenberg and cubic critical 
exponents are very small [2H], for example v differs by less than 0.1%, which 
is much smaller than the typical experimental error for Heisenberg systems 
[2]. Therefore, distinguishing the cubic and the Heisenberg universality class 
should be very hard in experiments and numerical Monte Carlo simulations. 

The stability properties of the 0(iV)-symmetric FP in the theory defined 
by Tic can be extended to any iV-component system in the presence of spin- 
4 terms. The point is that the cubic-symmetric perturbation J2i is a 
particular combination of the spin-4 operators Pf^ d and of the spin-0 term 
P 4 . Since the RG dimension of the spin-4 perturbation does not depend on 
the particular component, we can conclude that the 0(iV)-symmetric FP is 
unstable under any spin-4 quartic perturbation for iV > 3. 

4 Randomly dilute spin models 

The critical behavior of systems with quenched disorder is of considerable 
theoretical and experimental interest. A typical example is obtained by mix- 
ing a magnetic material with short-range interactions with a nonmagnetic 
material. These physical system can be described by the three-dimensional 
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Figure 1: Scketch of the RG flow in the coupling plane (u,v) for (a) Ising 
(M = 1) and (b) M-component (M > 1) randomly dilute systems. 



randomly dilute spin model 

= J Pi Pi s i ' s i> ( 23 ) 

<ij> 

where the sum is extended over all nearest-neighbor sites, are M-component 
spin variables, and pi are uncorrelated quenched random variables, which are 
equal to one with probability p (the spin concentration) and zero with prob- 
ability 1 — p (the impurity concentration). See, e.g., Refs. jU EUJ 12 HH] for 
reviews discussing these systems. For sufficiently low spin dilution 1 — p, 
i.e. above the percolation threshold of the spin concentration, the system 
described by the Hamiltonian 7i p undergoes a second-order phase transition 
atT c (p)<T c (p=l). 

The critical behavior of randomly dilute spin systems is rather well es- 
tablished. Randomly dilute Ising systems, of which the random Ising model 
(RIM) (J23|) is an example, belong to a new universality class which is distinct 
from the Ising universality class describing the critical behavior of the pure 
system. This has been clearly observed in experiments on dilute uniaxial an- 
tiferromagnets such as Fea.Zn.1_a.F2 and Mn^Zni-a^ materials. The critical 
exponents turn out to be independent from the impurity concentration and 
definitely different from those of the pure Ising universality class, see, e.g., 
their estimates reported in Refs. |3U1 l_l EH]. In these antiferromagnetic sys- 
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terns the presence of a uniform magnetic field causes a crossover to a different 
random- field critical behavior 32.. These experimental result are confirmed 
by renormalization-group studies, showing that the pure Ising FP is unstable 
with respect to the addition of impurities and that the critical behavior is 
controlled by a new stable RIM FP. 

Unlike Ising-like systems, multicomponent 0(M)-symmetric spin systems 
do not change their asymptotic critical behavior in the presence of random 
impurities. This is predicted by the Harris criterion |33] which states that the 
addition of impurities to a system which undergoes a continuous transition 
does not change its critical behavior if the specific-heat critical exponent 
a of the pure system is negative. For M > 1 the specific-heat exponent 
cum of the pure system is negative, for example 16, 17j a 2 ~ —0.014 and 
a% ~ —0.13, respectively for the XY and Heisenberg universality classes, and 
thus disorder is irrelevant. From the RG point of view, the Wilson-Fisher 
FP of the pure O(M) theory is stable under the perturbation induced by the 
dilution. The presence of impurities affects only the approach to the critical 
regime, leaving the asymptotic behavior unchanged. This is confirmed by 
several experiments that investigated the effect of disorder on the A-transition 
of 4 He (it belongs to the XY universality class) jHH ESI 12] and on Heisenberg 
magnets |H1 E3 EE] • 

The FT approach is based on an effective translation-invariant LGW 
Hamiltonian that is obtained by using the replica method jSHl HOI EU EE] , 
i.e. 

Hmn = \d d xW\ [(d^ ai f + ri&] + £ ± (u + v^j) ^% \ , (24) 

ia ijab ' ) 

where a, b = 1, ...M and i,j = 1, ...N. In the limit N — > the Hamiltonian 
Hmn with u < and v o > is expected to describe the critical properties 
of dilute M-component spin systems. Thus, their critical behavior can be 
investigated by studying the RG flow of TCmn in the limit iV — > 0. For 
generic values of M and N, the Hamiltonian Hmn describes M coupled iV- 
vector models and it is usually called MN model 4j. Figure Q sketches the 
expected flow diagram in the quartic-coupling plane, for Ising (M = 1) and 
multicomponent (M > 1) systems in the limit N — ■> 0. There are four FP's: 
the trivial Gaussian one, an 0(M)-symmetric FP, a self-avoiding walk (SAW) 
FP, and a mixed FP. The SAW FP is stable and corresponds to the (M x N)- 
vector theory for N — > 0; but it is not of interest for the critical behavior of 
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randomly dilute spin models, since it is located in the region u > 0, while 
the region relevant for quenched disordered systems corresponds to negative 
values of the coupling u. The stability of the other FP's depends on the value 
of M. Nonperturbative arguments OH] show that the stability of the O(M) 
FP is related to the specific-heat critical exponent of the 0(M)-symmetric 
theory. Indeed, TCmn at the 0(M)-symmetric FP can be interpreted as the 
Hamiltonian of N M- vector systems coupled by the 0(MiV)-symmetric term. 
Since this interaction is the sum of the products of the energy operators of 
the different M-vector models, the crossover exponent associated with the 
0(MiV)-symmetric quartic interaction is given by the specific- heat critical 
exponent % of the M-vector model, independently of N. This implies that 
for M = 1 (Ising-like systems) the pure Ising FP is unstable since <fi = a\ > 0, 
while for M > 1 the O(M) FP is stable given that dim < 0, in agreement 
with the Harris criterion. For M > 1 the mixed FP is in the region of positive 
u and is unstable 0J. Therefore, the RG flow of the M-component model 
with M > 1 is driven towards the pure O(M) FP. Thus, the asymptotic 
behavior remains unchanged when impurities are added. But their effect 
is not totally negligible because they give rise to slowly- decaying scaling 
corrections proportional to t£ with A r = — om, where a a/ is the specific- 
heat exponent of the pure 0(M)-symmetric theory. For Ising-like systems, 
the pure Ising FP is instead unstable, and the flow for negative values of 
the quartic coupling u leads to the stable mixed or random (RIM) FP which 
is located in the region of negative values of u. The above picture emerges 
clearly in the framework of the e expansion, although the RIM FP is of order 
\fe ^3] rather than e. 

The most precise FT results have been obtained in the framework of the 
FD expansion in powers of the zero-momentum quartic couplings related 
to u and Vq. In this scheme the theory is renormalized by introducing a 
set of zero-momentum conditions for the one-particle irreducible two-point 
and four-point correlation functions, which relate the zero-momentum mass 
scale m and the quartic couplings u and v to the corresponding Hamiltonian 
parameters r, uq, and Vq. This is just a straightforward extension of the 
method outlined in Sec. 12.11 to a theory with two quartic parameters. The 
critical behavior of dilute spin systems is determined by the stable FP of the 
theory, that is given by the common zero u*, v* of the /3-functions in the limit 
iV -> 



at \ 9U 

fJ u {u,v) = m— 



p v (u,v) = m— 

n , dq 



(25) 

u ,v 
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Figure 2: RG trajectories in the u, v plane for various values of s = uo/v < 
0. From Ref. 

whose stability matrix has positive eigenvalues (actually a positive real part 
is sufficient). Then, the critical exponents are obtained by evaluating appro- 
priate RG functions at u*,v*. Figs. 121 and 01 show typical RG trajectories in 
the zero-momentum quartic plane u, v in the relevant region corresponding 
to Mo < 0, for the RIM and the random Heisenberg model respectively 

The FD expansions of the (3 functions and of the critical exponents have 
been computed to six loops jlBJ I2H HZj • The results of the six-loop analyses 
are reported in Table El where they are compared with estimates obtained in 
Monte Carlo simulations of the RIM, see, e.g., Refs. jlHlESj, and in exper- 
iments on uniaxial antiferro magnets, such as Fe I Zni_ :r F2 and Mn x Zni_ x F2, 
see, e.g., Refs. j3Ul EI]- The agreement is satisfactory, although we note 
the slightly smaller value of the experimental estimate of a. We mention 
that several computations have also been done in the framework of the e 
expansion. The y/e expansion [H] turns out not to be effective for a quanti- 
tative study of the RIM (see, e.g., the analysis of the five-loop series done in 
Ref. jSDl)- The related minimal-subtraction renormalization scheme without 
e expansion $ : have been also considered. The three- and four-loop results 
turn out to be in reasonable agreement with the estimates obtained by other 
methods, but at five loops no random FP is found |5T]- Nonperturbative ap- 
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Figure 3: RG trajectories of the dilute Heisenberg model for various values 
of s = Uq/vq < 0. From Ref. |3E] . 



Table 3: Critical exponents for the RIM universality class. 
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V 


a 


P 


six-loop FD 


m 


1.330(17) 


0.678(10) 


-0.034(30) 


0.349(5) 


Monte Carlo 


HHJ 


1.342(6) 


0.683(3) 


-0.049(9) 


0.3534(15) 


Fe x Zni_ a; F2 


HD 


1.31(3) 


0.69(1) 


-0.10(2) 


0.359(9) 



proaches based on scaling-field expansions and on continuous RG equations 
have also been exploited [52J EHj > providing substantially consistent results. 
A more complete list of theoretical and experimental results can be found in 
Refs. EDI ED. 

Using the FT approach, one can also compute the critical exponent 
describing the crossover from random-dilution to random-field critical be- 
havior in Ising systems, and in particular the crossover observed in dilute 
anisotropic antiferromagnets when an external magnetic field is applied [3U] . 
The crossover exponent <fi is related to the RG dimensions of the quadratic 
operator (i ^ j) in the limit iV — > [Slj. Six-loop computations [55] 

provide the estimate <fi = 1.42(2), which is in good agreement with the avail- 
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able experimental estimates, for example <fi = 1.42(3) for Fe a; Zn 1 _ ; j,F2 |3*Uj . 

Finally, we mention that a perturbative six-loop analysis of the combined 
effect of impurities and magnetic anisotropy has been reported in Ref. [2H|- 



5 Frustrated spin models with noncollinear 
order 

In physical magnets noncollinear order is due to frustration that may arise 
either because of the special geometry of the lattice, or from the competition 
of different kinds of interactions jSHj- Typical examples of systems of the 
first type are stacked triangular antiferromagnets (STA's), where magnetic 
ions are located at each site of a three-dimensional stacked triangular lattice. 
On the basis of the structure of the ground state, in an iV-component STA 
one expects a transition associated with the breakdown of the symmetry 
from O(N) in the HT phase to 0(iV — 2) in the LT phase. The nature of 
the transition, and in particular the existence of a new chiral universality 
class |57j . is still controversial. On this issue, there is still much debate, FT 
methods, Monte Carlo simulations, and experiments providing contradictory 
results in many cases (see, e.g., Ref. [2] for a recent review of results). Overall, 
experiments on STA's favor a continuous transition belonging to a new chiral 
universality class. 

The determination of an effective LGW Hamiltonian describing the crit- 
ical behavior leads to the 0(M)®0(iV)-symmetric theory 



n ch = J ^^E[(W 2 +r^]+i«(E0a)V^E[(^-^) 2 -^ 

a, b 



(26) 

where a , a = 1, . . . , M, are iV-component vectors. The case M = 2 with 
v > describes frustrated spin models with noncollinear order. Negative 
values of v correspond to magnets with sinusoidal spin structures. XY and 
Heisenberg systems correspond to iV = 2 and N = 3 respectively. Recently 
the Hamiltonian (J26)) has also been considered to discuss the phase diagram 
of Mott insulators |55j . See Refs. [2)112] for other applications. 

Six-loop calculations [SH] in the framework of the d = 3 FD expansion 
provide a rather robust evidence for the existence of a new stable FP in the 
XY and Heisenberg cases corresponding to the conjectured chiral universality 
class, contradicting earlier studies based on much shorter (three-loop) series 
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60j. It has also been argued jnij that the stable chiral FP is actually a focus, 
due to the fact that the eigenvalues of its stability matrix turn out to have 
a nonzero imaginary part. The new chiral FP's found for N = 2,3 should 
describe the apparently continuous transitions observed in STA's. The FT 
estimates of the critical exponents are in satisfactory agreement with the 
experimental results, including the chiral crossover exponent related to the 
chiral degrees of freedom 

On the other hand, other FT studies, see, e.g., Refs. [HHJini|, based on 
approximate solutions of continuous RG equations, do not find a stable FP, 
thus favoring a weak first-order transition. Monte Carlo simulations have 
not been conclusive in setting the question, see, e.g., Refs. [123 EH1 EZ] • Since 
all the above approaches rely on different approximations and assumptions, 
their comparison and consistency is essential before considering the issue 
substantially understood. 

Finally, we mention that the six-loop FT analysis of the critical behavior 
of systems with nonplanar ordering, i.e. with M > 3, has not find any 
evidence of the presence of a stable FP (HE] , suggesting a first-order transition. 
A high-order FT study of two-dimensional systems have been reported in 
Ref. PI- 

6 Competition of two distinct order parame- 
ters 

The competition of distinct types of ordering gives rise to multicritical be- 
havior. More specifically, a multicritical point (MCP) is observed at the 
intersection of two critical lines characterized by different order parameters. 
MCP's arise in several physical contexts. The phase diagram of anisotropic 
antiferromagnets in a uniform magnetic field H\\ parallel to the anisotropy 
axis presents two critical lines in the temperature- Hn plane, belonging to the 
XY and Ising universality classes, that meet at a MCP [TUl [3j . MCP's are 
also expected in the temperature-doping phase diagram of high-T c supercon- 
ductors. Within the SO(5) theory [7TJ[72] of high-T c superconductivity, it has 
been speculated that the antiferromagnetic and superconducting transition 
lines meet at a MCP in the temperature-doping phase diagram. 

Different phase diagrams have been observed close to a MCP. If the tran- 
sition at the MCP is continuous, one may observe either a bicritical or a 
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tetracritical behavior. A bicritical behavior is characterized by the presence 
of a first-order line that starts at the MCP and separates the two different 
ordered low-temperature phases, see Fig. HJ In the tetracritical case, there 
exists a mixed low-temperature phase in which both types of ordering coexist 
and which is bounded by two critical lines meeting at the MCP, see Fig. |S] 
It is also possible that the transition at the MCP is of first order. A possible 
phase diagram is sketched in Fig. El In this case the two first-order lines, 
which start at the MCP and separate the disordered phase from the ordered 
phases, end in tricritical points and then continue as critical lines. 



6.1 Multicritical behavior in 0(ni)©0(n2) theories 

The multicritical behavior arising from the competition of two types of or- 
dering characterized by O(n) symmetries is determined by the RG flow of 
the most general 0(ni)©0(n2)-symmetric LGW Hamiltonian involving two 
fields <pi and <p 2 with rii and n 2 components respectively, i.e. |3] 

n mc = ( d^lUdM 2 + Hd^) 2 + \r l( f>l + \r 2 <j>\ (27) 



+u 1 {4>\) 2 + u 2 {4>\) 2 + w4>\4>l 

The critical behavior at the MCP is determined by the stable FP when both 
ri and r 2 are tuned to their critical value. An interesting possibility is that 
the stable FP has O(iV) symmetry, N = rii + n 2 , so that the symmetry gets 
effectively enlarged approaching the MCP. As we shall see, this is realized 
only in the case N = 2, i.e. when two Ising lines meet. 

The phase diagram of the model with Hamiltonian (J27|) has been inves- 
tigated within the mean- field approximation in Ref. [73]. If the transition 
at the MCP is continuous, one may observe either a bicritical or a tetracrit- 
ical behavior. But it is also possible that the transition at the MCP is of 
first order. A low-order calculation in the framework of the e expansion [3] 
shows that the isotropic 0(iV)-symmetric FP (iV = n\ + n 2 ) is stable for 
N < N c = 4 — 2e + 0(e 2 ). With increasing N, a new FP named biconal FP 
(BFP), which has only 0(n!)©0(r22) symmetry, becomes stable. Finally, for 
large N, the decoupled FP (DFP) is the stable FP. In this case, the two or- 
der parameters are effectively uncoupled at the MCP. The extension of these 
0(e) results to three dimensions suggests that for n\ — 1 and n 2 = 2, the 
case relevant for anisotropic antiferromagnets, the MCP belongs to the 0(3) 
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0(1^) 



T 



Figure 4: Phase diagram in the plane T-g presenting a bicritical point. Here, 
T is the temperature and g a second relevant parameter. The thick line ( "flop 
line" ) represents a first-order transition. 

universality class, while for rii = 2 and n 2 = 3, of relevance for the SO (5) 
theory of high-T c superconductivity, the stable FP is the BFP. 

The 0(e) computations provide useful indications on the RG flow in three 
dimensions, but a controlled extrapolation to e = 1 requires much longer se- 
ries and an accurate resummation exploiting their Borel summability. This 
has been recently achieved, extending the e expansion to 0(e 5 ) [231- A ro- 
bust picture of the RG flow predicted by the 0(rzi)©0(n2)-symmetric LGW 
theory has been achieved by supplementing the analysis of the e series with 
the results for the stability of the 0(N) FP (cf. Sec. OJ), which were also ob- 
tained by analyzing six-loop FD series, and with nonperturbative arguments 
allowing to establish the stability of the DFP [7U [75] . 

In order to establish the stability properties of the 0(N) symmetric 
FP, we first note that the Hamiltonian (|2*7)l contains combinations of spin- 
2 and spin-4 polynomial operators that are invariant under the symmetry 
0(ni)©0(n 2 ). Explicitly, they are given by 



V 2)0 = $ 2 , 



^2 = f>r 2 =^-> 2 > 
a=i iv 



(28) 



^4,0 = ($ 2 ) 2 , 



^4,2 = $ 2 P 2 ,2, 
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Figure 5: Phase diagram with a tetracritical point. 



4 " i3 Jki " " 2 ^ + (Af + 2)(JV + 4)' 

where $ is the A^-component field (0i,0 2 )- The RG dimension 1/2,2 of ^2,2, 
and therefore of the operator P 2 a 2 ( see Table EJ), provides the crossover expo- 
nent = t/2,2^ at the MCP. Concerning the fourth-order terms, the results 
for the spin-4 RG dimension 2/44 at the O(N) FP (see again Table |2J imply 
that the O(N) FP is stable only for N = 2, i.e. when two Ising-like critical 
lines meet. It is unstable in all cases with N > 3. This implies that for 
N > 3 the enlargement of the symmetry 0(ni)©0(n2) to 0(N) does not 
occur, unless an additional parameter is tuned beside those associated with 
the quadratic perturbations. 

The analysis of the five-loop e series shows that for N = 3, i.e. for ri\ = 1 
and ri2 = 2, the critical behavior at the MCP is described by the BFP, 
whose critical exponents turn out to be very close to those of the Heisenberg 
universality class. For N > 4 and for any ni,n 2 the DFP is stable, implying 
a tetracritical behavior. This can also be inferred by using nonperturbative 
arguments [73J [7S] that allow to determine the relevant stability eigenvalue 
from the critical exponents of the O(rij) universality classes. At the DFP, the 
dominant perturbation is the quartic coupling term w<p\(f)\ that scales as the 
product of two energy-like operators, which have RG dimensions (1 — a^jvi 
where «j and z/j are the critical exponents of the O(rij) universality classes. 
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Figure 6: Phase diagram with a first-order MCP. The thick lines represent 
first-order transitions. 



Therefore, the RG dimension related to the ^-perturbation is given by 

«i aco 1 1 , 
y w = ^ + ^ = - + --d. (29) 

The stability properties of the DFP can be established by determining the 
sign of y w . Using the estimates of the critical exponents of the three-dimension- 
al O(rij) universality classes, see Sec. EJ y w turns out to be negative for 
N = rt\ +7i2 > 4, and positive for = 2, 3. Therefore, the DFP turns out to 
be stable for N > A. Finally, when the initial parameters of the Hamiltonian 
are not in the attraction domain of the stable FP, the transition between the 
disordered and ordered phases should be of first order in the neighborhood 
of the MCP. In this possible phase diagram is given in Fig. |Hl where 

all transition lines are first-order close to the MCP. 



6.2 Anisotropic antiferromagnets in a uniform mag- 
netic field 

Anisotropic antiferromagnets in a uniform magnetic field Hn parallel to the 
anisotropy axis present a MCP in the T — Hn phase diagram, where two 
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critical lines belonging to the XY and Ising universality classes meet [5J. In 
this case the high-order FT analyses predict a multicritical behavior described 
by the BFP, in contrast with earlier low-order e-expansion calculations [3] 
that predicted the enlargement of the symmetry to 0(3). The mean-field 
approximation assigns a tetracritical behavior to the MCP [3J, but a more 
rigorous characterization, that requires the computation of the corresponding 
scaling free energy, is needed to draw a definite conclusion. Experimentally, 
the MCP appears to be bicritical, see, e.g., Refs. jTHl [77] . A quantitative 
analysis of the biconal FP shows that its critical exponents are very close to 
the Heisenberg ones. For instance, the correlation-length exponent v differs 
by less than 0.001 in the two cases. Thus, it should be very hard to distinguish 
the biconal from the 0(3) critical behavior in experiments or numerical works 
based on Monte Carlo simulations. The crossover exponent describing the 
crossover from the 0(3) critical behavior is very small, i.e., = y^v ~ 
0.01, see TableEl so that systems with a small effective breaking of the 0(3) 
symmetry show a very slow crossover towards the biconal critical behavior 
or, if the system is outside the attraction domain of the BFP, towards a 
first-order transition; thus, they may show the eventual asymptotic behavior 
only for very small values of the reduced temperature. 

6.3 Nature of the multicritical point within the SO(5) 
theory of high-T c superconductors. 

High-T c superconductors are other interesting physical systems in which 
MCP's may arise from the competition of different order parameters. At 
low temperatures these materials exhibit superconductivity and antiferro- 
magnetism depending on doping. The SO (5) theory [7TJ [72] attempts to 
provide a unified description of these two phenomena, by introducing a three- 
component antiferromagnetic order parameter and a <i-wave superconducting 
order parameter with U(l) symmetry, with an approximate 0(5) symmetry. 
This theory predicts a MCP arising from the competition of these two order 
parameters when the corresponding critical lines meet in the temperature- 
doping phase diagram. Neglecting the fluctuations of the magnetic field and 
the quenched randomness introduced by doping, see, e.g., Ref. [73] for a 
critical discussion of this point, one may consider the 0(3)©0(2)-symmetric 
LGW Hamiltonian to infer the critical behavior at the MCP, see, e.g., Refs. 
UHl E3 E01 IHH IH21 E3 In particular, the analysis of Ref. [HI], which uses the 



22 



projected SO(5) model as a starting point, shows that one can use Eq. ()27|) 
as an effective Hamiltonian. 3 

Different scenarios have been proposed for the critical behavior at the 
MCP. In Refs. [7111111113, it was speculated that the MCP is a bicriti- 
cal point where the 0(5) symmetry is asymptotically realized. This picture 
has been apparently supported by Monte simulations for a five-component 
0(3)©0(2)-symmetric spin model [7jl|, and by a quantum Monte Carlo study 
of the quantum projected SO (5) model in three dimensions |8%| . These 
numerical studies show that, within the parameter ranges considered, the 
scaling behavior at the MCP is consistent with an 0(5)-symmetric critical 
behavior. On the basis of the 0(e) results of Refs. Refs. jH01IH21 predicted 
a tetracritical behavior governed by the BFP. However, since the BFP was 
expected to be close to the 0(5) FP, it was suggested that at the MCP the 
critical exponents were in any case close to the 0(5) ones. These hypothe- 
ses are contradicted by the high-order FT results that we have discussed 
in Sec. 16.11 They indeed show that the DFP is the only stable FP. This 
implies that the asymptotic approach to the MCP is characterized by a de- 
coupled tetracritical behavior or by a first-order transition for systems that 
are outside the attraction domain of the DFP. The 0(5) symmetry can be 
asymptotically realized only by tuning a further relevant parameter, beside 
the double tuning required to approach the MCP. The 0(5) FP is unstable 
with a crossover exponent (f)^ = y^v f& 0.15, see Table El which, although 
rather small, is nonetheless sufficiently large not to exclude the possibility of 
observing the RG flow towards the eventual asymptotic behavior for reason- 
able values of the reduced temperature. Evidence in favor of a tetracritical 
behavior has been recently provided by a number of experiments, see, e.g., 
Refs. jHHl EH IH7J EH1 EHl EOj , which seem to show the existence of a coex- 
istence region of the antiferromagnetic and superconductivity phases. The 
possible coexistence of the two phases has been discussed in Refs. jHSl EH E2] • 
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